Geometric potential and transport in photonic topological crystals 

Alexander Szameit^*, Felix Dreisow^, Matthias Heinrich^, Robert 
KeiP, Stefan Nolte^, Andreas Tiinnermann^, Stefano Longhi^ 
^Physics Department and Solid State Institute, Technion, 32000 Haifa, Israel 
"^Institute of Applied Physics, Friedrich-Schiller-Universitdt, 
Max-Wien-Platz 1, 07743 Jena, Germany and 
^ Dipartimento di Fisica, Politecnico di Milano, 
Piazza Leonardo da Vinci 32, 20133 Milan, Italy 

Abstract 

We report on the experimental reahzation of an optical analogue of a quantum geometric poten- 
tial for light wave packets constrained on thin dielectric guiding layers fabricated in silica by the 
femtosecond laser writing technology. We further demonstrate the optical version of a topological 
crystal, with the observation of Bloch oscillations and Zener tunneling of purely geometric nature. 
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Studying the quantum mechanics of a particle constrained on a curved space jl-^ has 
been a matter of great controversies for more than 50 years. For a long time, it has been 
known that geometry influences in different ways the motion of classical and quantum parti- 
cles confined on a curved surface or on a line. When surface curvatures become comparable 
to the de Broglie wavelength, geometrical effects come into play by the internal metric of the 
constraint manifold as well as by the external metric of the embedding space. As opposed 
to a classical particle, a quantum particle retains some knowledge of the surrounding three- 
dimensional space and, in spite of the absence of interactions, it experiences an effective 
frictional potential of geometric nature {2, 3|. Such a geometric potential is of major impor- 
tance in the understanding and control of the physical properties of novel low-dimensional 
functional materials, like curved carbon nanotubes and DNA wires j4]. For instance, an 
electron confined to a periodically curved surface senses a periodic frictional potential which 
acts as a topological crystal The definition of a geometric potential in quantum and 
condensed matter physics, however, is troubled by the realization of a proper squeezing 

na nn 

procedure p, |7|, which is needed to avoid operator-ordering ambiguities [8|, |9[. A rather 
accepted approach that avoids operator-ordering ambiguities is the confining-potential ap- 
proach, originally proposed by Jensen, Koppe and da Costa (JKC) 2, sl, in which a strong 
force acting normal to the surface provides the appropriate confinement. Quantum excita- 
tion energies in the normal direction are raised far beyond those in the tangential direction. 
Hence, the particle motion normal to the surface can be safely neglected, which leads to an 
effective Hamiltonian for propagation along the curved surface. Though the JKC approach 
has gained a broad consensus especially in the theoretical condensed matter physics com- 
munity [lOHlJ], it was noticed that such an ideal squeezing procedure may correspond to 
unrealistic restrictions or cannot be unique 

a a. 

Clear evidences of geometric potentials 
in low- dimensional curved nanosystems are still lacking and are likely to be a controversial 
matter. On the other hand, optics has offered in recent years a fascinating laboratory tool to 
investigate classical analogues of otherwise inaccessible quantum-mechanical and relativistic 
effects (see, e.g., 15l-ll9|). In this Letter we report on the first experimental realization of an 
optical analogue of a quantum geometric potential for light wave packets constrained on thin 
dielectric guiding layers fabricated in silica by the femtosecond laser writing technology |20 |. 



We further demonstrate the optical version of a topological crystal [5| , with the observation 
of Bloch oscillations and Zener tunneling of purely geometric nature |2l[ |. 
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Let us consider wave propagation at wavelength A in a thin and weakly guiding dielectric 
layer, which is assumed to be invariant along the z direction and arbitrarily curved in the 
transverse (x, y) plane along a curve 7 (see Fig. la). In the scalar approximation, the electric 
field E{x,y,z) satisfies the Helmholtz equation 

V',E + d',E+{n/XrE = 0, (1) 

where = 9^ + dy is the transverse Laplacian, A = A/ (27r) is the reduced wavelength, and 
n = n{x, y) is the refractive index profile of the optical structure, which is assumed to be 
independent of z and weakly deviating from the refractive index Ug of the substrate. The 
analogy between spatial propagation of light waves along the z direction and the temporal 
dynamics of a two-dimensional quantum particle constrained on the curved line 7, discussed 
by da Costa [3], is at best captured by the introduction of the local curvilinear coordinates 
((J, 77) of Fig. la. The volume element in the curvilinear system is given by dV = xdcrdrjdz, 
where ^7) = I + rj/R and R = R{cr) is the local radius of curvature of 7. To study the 
behavior of light waves near the guiding layer (77 — )■ 0), following 3| it is worth introducing 
the new wave field S = E^fx- After writing in Eq.(l) the transverse Laplacian Vj in 
curvilinear coordinates and with the substitution E = £ / ^Jxi the evolution equation of the 
field £ in the limit — (i.e., close to the line 7) reads 

{dl + dl + dl)£ + {IIAR')£ + {nlXf£ = 0. (2) 

In the weak guiding limit |72(cr, r]) — ns\ ^ and assuming \R\ ^ A, the paraxial approxi- 
mation can be introduced in the usual way by setting £{a, t], z) = ipio', V: ^) exp(ms-2/A) and 
neglecting dlip as compared to {l/\)dz'ip. This yields the optical Schrodinger equation 

= -7^ + l£ V + [y^ia, V) + V,ia)] ^, (3) 



dz 2ns xda"^ drj 
where Vc{(J, rj) = [nl — n'^{a, ri)]/{2ns) ~ — n(cr, rj) is the confining potential and 



is the so-called geometric potential [3]. The correspondence of Eq.(3) to the quantum- 
mechanical Schrodinger equation of a two-dimensional particle constrained on the curve 7, 
discussed by da Costa , is formally obtained after replacing the photon wavelength A with 
the Planck constant h, the refractive index Ug with the particle mass m and the spatial 



coordinate z with time t. Hence, in optics the quantum-mechanical evolution in time of a 
two-dimensional wave function is mapped onto the propagation of an optical wave packet 
along the spatial 2;-direction. The ^-independence of the path 7 reflects the circumstance 
that the constraint is time-independent. Similarly to the quantum mechanical problem jsl, 
in the optical Schrodinger equation (3) the confining potential Vc, squeezing the wave packet 
around 7, originates from the physical change of the refractive index in the guiding layer 
(see Fig. lb); on the contrary, the geometric potential Vg is of geometric nature and arises 
from the diffraction operator — (A^/2?2s)Vj (the analogue of the kinetic energy operator) 
in curvilinear coordinates. In the curved reference frame (ex, 77), this potential acts on light 
waves like a fictitious refractive index change. 

In the JKC approach, the confining potential Vc is taken to be independent of a, a condition 
which ensures that the constraint is frictionless in the classical (geometric-optic) limit. In 
this case, the motion in the normal t] direction can be exactly separated from the dynamics by 
letting ip{a, i], z) = F{a, z)g{ri) exp{—iEoz/X), where g{ri) is the ground-state wave function 
of the confining potential Vc and Eq its corresponding energy [i.e., —{X^/2ns){cPg/drf) + 
Vc{ri)g{ri) = Eog^r])]. One then obtaines an effective equation for the particle motion along 
the a coordinate 

Therefore, the geometric potential V„ acts as an effective frictional potential for the motion 



on the curved manifold. Unfortunately, t 
potentials have this exceptional property 



lere is no compelling reason that actual confining 
6|. Moreover, it was shown that one can construct 



ad-hoc confining potentials leading to an arbitrary correction to the geometric potential 
that obey the frictionless condition of the constraint in the classical limit A — )■ [7|. Non- 
ideal squeezing potentials will generally replace the geometric potential Vg in Eq.(5) by an 
effective frictional potential K//- If we allow the confining potential V^ to slowly vary with 
cr, the tangential motion can still be approximately separated from the normal one by using 
a multiple scale asymptotic analysis (as in Q]). In this way, the tangential motion turns 
out to be governed again by the reduced Schrodinger equation (5), but with the geometric 
potential Vg replaced by V^fj = Vg + AV^, where 

^""'^"^ = ^Eia) - ^ / A«0. (6) 
In the above equation, g{a, rj) is the local ground-state wave function of the confining poten- 



tial Vc, normalized such that J dr]g^ = 1, and Eq + AE{(t) is the corresponding cr-dependent 
energy, i.e. —{X^/2ns)d'^g + VcQ = {Eq + A.E)g. The JKC geometric potential is attained 
whenever AVg{a) vanishes, a condition which does not necessarily imply the cr-invariance of 

In order to demonstrate the features of an ideal JKC squeezing procedure and the dis- 
crepancies of a non-ideal squeezing, we realized an optical analog of a topological crystal sl, 
in which a two-dimensional wave function is squeezed onto an undulating curve 7. To this 
aim, a sinusoidally undulated slab waveguide (Fig. [1};), defined by y{x) = A sin(27rx/A), 
was fabricated using the laser direct-writing technology 201]. The use of this technique com- 
monly yields a layer which exhibits constant thickness in the vertical y-direction rather then 
in the normal r^-direction. For a homogeneous refractive index change as shown in Fig. [T]i, 
this implies that the correction to the geometric potential does not vanish, i.e. AVg 7^ 0. 
However, using the laser direct-writing approach, the refractive index variations along 7 
can be tailored with great accuracy. In particular, it is possible to modulate the strength 
of n{a,ri) along a such that AVg gets negligible as compared to Vg. The refractive index 
distribution that closely approximates the JKC condition AVg = is shown in Fig. [1^, and 
corresponds to the highest index change at the inflection points of the undulation, where 
R = 00. 

To investigate wave packet dynamics in the topological crystal, we employ a fluorescence 
microscopy technique 22] . Light around A = 633 nm excites color centers, which are formed 
during the waveguide fabrication process. The resulting fluorescence is proportional to the 
intensity of the propagating light, and can be observed from above the sample, as sketched 
in Fig. [21 due to the isotropic emission. When launching a beam into the two samples 
with the refractive index profiles representing the non-ideal (Fig. [T]i) and the ideal (JKC) 
case (Fig. [1^), one obtains the intensity distributions shown in Fig. [3^,c, respectively. The 
observations are confirmed by numerical simulations of the full wave equation (see Fig. [3b, d). 
Although in both cases the initial wave packet spreads, the spreading rate is larger for the 
JKC potential. Most importantly, the spreading pattern observed in Fig. 3c reproduces with 
excellent accuracy the diffraction pattern in a one-dimensional crystal with the potential 
defined by Eq.(4). Hence, our experimental technique enables to squeeze a wave packet 
mimicking the ideal JKC method. It is important to point out that the equivalent lattice 
potential Vg caused by the curvature of the JKC guide has a geometric origin and cannot be 
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naively exp lained by effective-index or variational methods generally adopted in guided-wave 



optics 



23| . As the effective index method well explains the formation of periodic optical 



potentials in common waveguide array settings (such as those investigated in 2^), it fails to 
explain the onset of the geometric potential for the JKC undulating slab structure of Fig.le. 
In fact, if one calculates from Eq.(l) the effective index along the vertical y direction using 
a separation variable method [23|, the resulting effective index ne{x) [25| turns out to be 
largest in the regions with vanishing curvature {R = oo), where the local refractive index 
is highest (see Fig.le). This would erroneously predict light confinement in regions around 
R = oo, rather than where \R\ is minimal [see Eq.(4)]. 

Using the confirmed JKC potential, we finally investigated the transport properties in 
the topological crystal by the application of a direct current (dc) force. Similarly to ordinary 
waveguide arrays, the dc force is expected to inhibit wave packet spreading and to cause 
an oscillatory motion of the wave packet [Bloch oscillations (BOs)] via the formation of a 
Wannier-Stark ladder spectrum [2J], with the occurrence of Zener tunneling (ZT) at high dc 
forcing 26 1. According to a fictitious dc force is realized by fabricating the undulated 
slab slightly curved in the ^-direction, which approximately preserves the z-invariance of the 
path 7, but imprints a transverse dc force in the x-direction. Our structure thus realizes a test 
bed for the observation of so-called topological BOs predicted in 2l|, where both effective 
lattice potential and dc force arise from the geometric deformation of a slab waveguide. The 
experimental observation of BOs in topological photonic crystals, corresponding to a broad 
input beam excitation, is shown in Fig. Hti, together with the numerical prediction (Fig. 
Hb). In order to visualize ZT, we doubled the transverse force by halving the longitudinal 
radius of curvature of the undulated slab waveguide. The results corresponding to broad 
beam excitation at the input plane are shown in Figs. |lt,d. Note that, as compared to Figs. 
UK andHb, ZT is now clearly visible because a considerable fraction of the wave packet does 
not experience Bragg reflection, but tunnels into higher propagation bands. 

In conclusion, we experimentally observed a geometric potential for optical wave packets 
constrained on curved surfaces and showed the impact of the squeezing procedure on the 
resulting frictional potential. In particular, an optical version of a topological crystal sl 
has been realized, in which the band structure of the crystal is determined by the geometric 
potential according to the JKC theory [3]. Our experiments shed new light onto the old and 
rather controversial problem of wave mechanics of quantum particles constrained on a curved 
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space, and open the possibility to explore the interplay between topology and transport 
in low-dimensional curved structures. Our findings could also pave the way towards the 
investigation of photonic materials with topologically-controUed diffractive and refractive 
properties. 
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FIG. 1: (a) Sketch of a two-dimensional cylindrical (i.e. z- invariant) curved surface with the 
curvilinear coordinates {a, ij). (b) The forced path of a wave packet in a constraining potential, (c) 
Microscope image of the front facet of the undulating guiding layer that realizes a one-dimensional 
topological crystal, (d) Homogeneous refractive index change yielding non-ideal squeezing, (e) 
Corrected refractive index distribution, resulting in ideal JKC conditions. 
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Experiment Simulation 




FIG. 3: (a) Measurement of the light propagation in a non-ideal squeezing potential, and (b) 
corresponding numerical confirmation, (c) Measurement of the light propagation under ideal JKC 
conditions, and (d) corresponding numerical confirmation. 
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Exp&rinneflt Simulation 




FIG. 4: (a) Dynamics of light in an optical topological crystal for broad beam excitation, when 
a transverse dc force is applied, (b) Numerical confirmation, (c) Zener tunneling for broad beam 
excitation in a high potential gradient, (d) Numerical confirmation. Note that the transverse dc 
force is applied by slightly curving the guiding layer in the z-direction. For reasons of visibility, this 
curvature has been numerically removed from the experimental images. Instead, the transverse 
gradient is sketched by tilting the images. 
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